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CHAPTER  I 


INTROnUCTION 

In  applying  the  Geometrical  Theory  of  Diffraction  (GTD)  to  antenna 
radiation  and  coupling  problems  involving  curved  surfaces.  A  major  task 
is  to  determine  the  geodesic  paths  on  the  curved  surface.  For  airborne 
antennas  mounted  on  the  fuselage  of  an  aircraft,  the  fuselage  is 
generally  modeled  as  a  cylinder  or  a  prolate  spheroid  [1,2]  in  the  GTD 
analysis.  However  to  better  approximate  general  fuselage  shapes,  an 
ellipsoidal  model  of  the  fuselage  is  needed.  Among  the  solutions  for 
obtaining  the  geodesic  paths  of  the  three  models,  i.e.,  the  cylindrical, 
the  prolate  spheroidal  and  the  ellipsoidal  models,  the  ellipsoid  case  is 
the  most  involved  and  complex  one.  This  is  to  be  expected  because  the 
equation  describing  an  ellipsoid  can  degenerate  into  that  of  a  cylinder 
or  a  prolate  spheroid  by  using  appropriate  parameters. 

According  to  the  generalized  Fermat's  principle,  a  ray  emanating 
from  a  source,  which  is  located  on  the  surface,  follows  a  geodesic  path 
on  the  surface  and  continually  sheds  energy  into  the  shadow  region. 

Such  a  creeping  wave  mechanism  is  illustrated  in  Figure  1,  from  which  it 
can  be  seen  that  a  ray  traverses  from  the  source  point  O'  to  the 
diffraction  point  Q,  and  then  propagates  along  the  geodesic  tangent  at  0 
toward  the  observation  point  Ps. 

As  the  energy  flows  around  the  surface,  it  is  continuously 
diffracted  along  the  geodesic  tangent  toward  the  field  point  such  that 
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Figure  1.  «ey  gatli  in  the  shado-  region. 
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the  significant  effect  of  the  surface  is  associated  with  a  region  around 
the  source.  In  fact,  for  an  ellipsoid,  the  significant  portion  of  the 
surface,  which  is  associated  with  the  dominant  energy,  may  look  as  shown 
in  Figure  2. 

Note  that  this  region  can  be  specified  by  following  the  various 
geodesic  paths  until  the  radiation  level  along  a  given  path  becomes 
insignificant,  i.e.,  more  than  40db  below  the  source  magnitude.  With 
this  in  mind,  it  is  clear  that  one  could  represent  the  ellipsoid  by  a 
structure  which  simulates  the  elliptic  cross  section  completely; 
however,  the  profile  could  be  approximated  by  a  simpler  shape  since  the 
significant  energy  region  does  not  cover  a  large  portion  of  the  profile 
shape.  An  elliptic  cone  model  is  employed  here  to  simulate  the 
ellipsoid  which  in  turn  can  be  used  to  model  a  fuselage.  This 
perturbation  model  is  illustrated  in  Figure  3  for  a  source  located  near 
one  end  of  the  ellipsoid.  Note  that  if  the  source  is  placed  at  the 
center  of  the  ellipsoid,  the  elliptic  cone  actually  becomes  a  right 
elliptic  cylinder. 

Since  the  elliptic  cone  is  a  developed  surface,  one  can  unfold  the 
elliptic  cone  such  that  a  planar  structure  results.  The  geodesics 
associated  with  the  elliptic  cone  are,  then,  straight  lines  on  this 
planar  structure.  In  order  to  allow  for  a  geodesic  solution  between  the 
simplicity  of  the  elliptic  cone  and  the  rigor  of  the  ellipsoid,  one  can 
perturb  the  elliptic  cone  by  bending  it  along  its  generator  as 
illustrated  in  Figure  3(b).  In  that  a  perturbation  technique  is  employ¬ 
ed,  the  geodesic  paths  for  the  elliptic  cone  are  simply  modified  such 
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ANTENNA 


REGION  OF 
SIGNIFICANT 
GEODESIC  PATHS 


Figure  2.  The  region  of  significant  energy  flow  from 
antenna  mounted  on  an  ellipsoid. 


that  the  solution  for  the  ellipsoid  is  basically  straight-forward  and  is 
the  basis  for  this  report.  It  is  obvious  that  one  cannot  use  this 
perturbation  technique  if  significant  energy  propagates  far  away  from 
the  source.  However,  as  mentioned  previously,  the  energy  which 
propagates  great  distances  along  the  ellipsoid  surface  becomes  insign¬ 
ificant  in  magnitude  such  that  one  need  not  solve  for  the  true  geodesic 
paths  outside  the  significant  region  shown  in  Figure  3(b).  The 
simplicity  of  these  perturbed  geodesic  paths  allows  one  to  very 
efficiently  determine  the  significant  ray  paths  on  the  ellipsoid. 


CHAPTER  II 

THEORETICAL  BACKGROUND 


A.  Introduction 

The  radiated  field  of  an  ellipsoid-mounted  antenna  is  analyzed 
using  the  Geometrical  Theory  of  Diffraction  (GTD).  The  surface  is 
assumed  to  be  perfectly  conducting,  and  the  surrounding  medium  is  free 
space.  An  exp  (jojt)  time  dependence  is  understood  and  suppressed  in  the 
following  formulations. 

Consider  an  infinitesimal,  magnetic  current  moment  dF„,(Q')  or  an 
electric  current  moment  dFg(Q')  located  on  a  perfectly  conducting  convex 
surface  as  shown  in  Figure  1;  the  sources 

dPn,(Q')  =  E(Q')xn'da' 

A 

and  c(Fe(Q')  =  I()t')dt'n' 
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pertain  to  the  aperture  and  monopole  type  excitations  with 


r(Q')  =  electric  field  at  Q', 

n'  =  outward  unit  surface  normal  at  Q', 

da'  =  area  element  at  Q', 

I(t')  =  electric  current  distribution  on  the  monopole,  and 
t'  =  distance  parameter  along  the  monopole. 

According  to  geometrical  optics,  the  space  surrounding  the  source  is 
divided  into  an  illuminated  and  shadow  region  by  a  plane  tangent  to  the 
surface  at  Q'.  This  plane  is  referred  to  as  a  shadow  boundary.  Since 
the  field  in  the  deep  lit  region  is  essentially  that  obtained  from 
geometrical  optics,  and  the  field  in  the  deep  shadow  region  is 
relatively  weak,  the  solution  for  the  transition  region  adjacent  to  the 
shadow  boundary  is  of  more  interest  and  discussed  below. 

B.  Shadow  Region 

The  creeping  wave  mechanism  in  the  shadow  region  is  illustrated  in 

Figure  4.  From  the  generalized  Fermat's  principle,  a  ray  emanating 

from  the  source  dPm(q')  at  Q'  traverses  a  geodesic  path  Q'Q  on  the 
e 

surface,  and  propagates  along  the  geodesic  tangent  at  Q  toward  the  field 

point  Pe.  The  field  dTm  at  Pc  can  be  expressed  in  terms  of  the  field  at 

e 

a  reference  point  Pg  by  [3] 

d  d  -jksg 

drm(Ps)'^rm(Po)  pip2 _  e  +0[m-2,m-3]  (i) 

6  0  d  d 

y  (pI+Sq)  (p2+So) 


where 


d  d 

PI  and  P2  are  the  principal  radii  of  curvatures  of  the  wavefront 
at  Pq; 

and  0[m*2,tn”3, . . .]  are  the  higher  order  terms. 

From  Figure  4,  it  is  seen  that  if  the  reference  point  Pq  is  moved  to  the 

d  d 

curved  surface  diffraction  point  Q,  then  pi+0,  pj+pc*  ^nd  Sg+s.  Since 
drm(Ps)  is  independent  of  the  reference  point  Pg,  it  follows  that 


li"'  /p?  dlm(Po)  =  rm(Q'.O)  ;  (2) 

^  ®  e 

Po-^ 

d 

Pl+0 


then. 


~jks 

dTm(Ps)~l]^(q-,  Q)  , Pc  e 
e  e  i' 's(pc+s) 


(3) 


Furthermore,  IiD(Q',Q)  can  be  related  to  the  source  strength  d^m  at  Q'  by 
e  ® 


ijn(Q'.Q)  =  dFm  (Q* ) -IrnCQ' ,Q) 

e  e  e 

3 

where  Tm(Q',Q)  is  given  by  [3] 
e 


(4) 


Tm(Q',Q)  =4^b'nT  (Q')H+t'bT  (q')S+b*bT  (Q')S+fnT  (0')H] 
4n  1  2  3  4 


>jkt 

e 


r  Pg(Q)i  1/6 


(5) 


Te(Q'.Q) 


[n‘nT  (Q')H+n'bT  (Q’)S] 
Tit  5  6 


-jkt 

e 


r  Pg(Q)i  1/6 


(6) 


Here  (t',n',b*)  and  (t,n,b)  are  the  tangent,  normal  and  binormal  unit 
vectors  to  the  surface  at  the  source  point  (Q')  and  diffraction  point 

AAA  AAA 

(Q),  respectively.  As  seen  from  Figure  4,  t  x  n  =  b  and  t'  x  n'  =  b'. 
The  quantities  Ti(Q'),  15(0')  are  the  torsion  factors  at  Q'  and  are 

given  in  Table  1.  Also, 


H  =  g(c) 

(7) 

(8) 

with 

g(?)  =  _1  /  dr  exp{-j rg)  and  (9) 

^  «exp(-j2ir/3)  w^(t) 
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TABLE  I 


(FOR  SHADOW  REGION) 


TYPE  OP  CONVEX 

SURFACE 

SLOT  o«  e*sc 

_ 

MONOPOLE  OR 

4  CASE 

_ 

SURFACE  RAY 

TORSfON 

TIO'I 

tURFACC  RADIUS 

OPCURVAYURC  IN 

t  DIRECTION 

SURFACE 

Diffracted 

RAY  CAUSTIC 

distance 

T,  (O'J 

T,  (O'l 

(0‘» 

TjlOl 

Tp'O’l 

SPHERE 

R 

■ 

0 

R 

R 

0 

0 

0 

0  t««(7) 

CIRCULAR 

CVtINOER 

■ 

R 

IHff 

R 

R 

•ift  la' 

• 

■ 

2«  a‘ 

2« 

•u*  a* 

ARBITRARY 

CONVEX 

SURPACf 

■ 

R 

TIO'l/»  lO'l 

R 

R 

no' \  filin') 

»■»««'  1  _2_ _ 

2  V  A, (0  1  R,ie  )/ 

WITH  R|IO’l  t  R|ia’l 

Note:  (1)  a'  is  defined  by  Tj.t'=coso'  where  rj  is  the  principal  direc¬ 


tion  unit  vector  associated  with  Ri(Q’). 


(2)  The  quantities  E  and  G  denote  two  of  the  three  coefficients 
E,  F,  G  that  appear  in  the  "first  fundamental  form"  of 
Differential  Geometry  [53. 
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g(^)  =  _L  J  dT  exp(-jTC)  (10) 

^  <-exp(-j2iT/3)  w^(t) 


which  are  known  as  the  acoustic  hard  and  soft  Fock  functions.  The  Fock 


type  Airy  function  is  given  by 


00 

W2(t)  =  _1_  /  dt*exp  (Tt-t^/3) 

^  «exp(j2iT/3) 


(11) 


and  w^(t)  is  the  derivative  of  W2(t)  with 
meter  ?  for  the  shadow  region  is  given  by 


respect  to  t. 
[3] 


The  Fock  para 


C  =  /  dt  ■ 


m(t') 

Pg(t  ) 


with 


m(t' ) 


kp  (f) 

g 

2 


1/3 


(12) 


(13) 


Here  Pg(t')  is  the  surface  radius  of  curvature  along  the  ray  path  at  t'. 
The  width  of  the  surface  ray  tube  at  Q,  dn(Q),  is  given  by 


dn(Q)  =  p  d.  .  (14) 

c  ^ 

Note  the  parameters  Zq  and  t  are  defined  as  the  free  space  wave  impedance 
and  geodesic  arc  length  from  O'  to  Q,  respectively. 
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Combining  Equations  (3)-(14),  the  n  and  b  directed  components  of 

drm(P  )  are  given  by  [3] 
e  s 


(a)  d'Pm(Q')  case: 


dE"(P  )  =  ^(dPm-b')He-jkt[3LS^]  /  o 

^  ^  pg{0)  y  3^ 

e"J‘^^+0[m‘^] 


s(pc+s) 


(15) 


A  A 

dE^(p5)  =  ^  [(dPm-b)TQS+(dPro-f)S] 


Pg(Q) 


^  Vs(pc+s) 


1  e-j'‘^0[m-2.m-3] 


(16) 


(b)  dPe(Q')  case: 


dE"(Ps)  =  dP  fnMHP~j*^^rPg(Q'h 

®  ®  PglHF 

e-j'^^OCm-2] 


“31  Vs(pc+s) 


(17) 


dEj(Ps)  *  dP, 


.(q')T  Se~j'^^[M^.?,]  >/^  J  .1  ~ 

°  Pg(Q)  V 


e-^'^SoCm-2] 


(18) 


where  To*T(Q' )pg(Q' )  with  T(Q')  being  the  surface  torsion  at  the  source 
location  (refer  to  Table  1). 
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C.  Lit  Region 


From  geometrical  optics,  the  source  dFm{Q')  at  O'  excites  waves 

e 

which  propagate  along  straight  line  ray  paths  from  the  source  to  field 

point  in  the  lit  region.  As  shown  in  Figure  5,  the  field  drm(P  )  at 

e  L 

point  Pl  is  expressed  by 


-jkSo 

+0[m-2,m-3]  (19) 

I  I 


Since  Q'  is  the  only  caustic  of  the  incident  rays,  the  principal  radii 
of  curvature  J  and  J  associated  with  the  incident  wavefront  at  Pr,  are 

Pl  P2  ^ 

identical,  i.e.,  p^=p^=  p^ •  Furthermore,  drm(P  )  is  independent  of  the 


reference  point  Pg.  If  Pg  is  chosen  to  be  at  O',  it  follows  that 


l 

lim  pidrm(Po)  =  Tm  (20) 

e  e 

Po-Q' 

pi  ■*■0 
So>S 

Ji 

should  exist.  Thus,  Lm  can  be  related  to  dFm(Q')  by  [3] 

e  e 


i. 

Tm  =  dFm(Q' ) ’fm  (21 ) 

e  e  e 
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Fquations  (19)-(21)  are,  then,  combined  to  yield 

-  -  T 

dEg{P|_)~<lP^(Q‘  — +0[m-'^,m--5, ...] 

=  Jlii  [S‘nA+t ‘  bB+b'ftc+^ '  nO] 

4^  t 

"^e  =  [^nM+n'^N] 

4  if 


(22) 

(23) 

(24) 


with  A,  B,  C,  D,  M,  and  N  are  defined  in  Table  2.  Note  that  qE^CP^^) 
decoupled  into  n  and  b  components  as  follows: 


(a)  dPn,(Q')  case: 


dE;;{PJ  =  ^  [(dP„-6')(HV^^FcoS«')  *  (dP„-{')1„r-C0S8'] 


*  o[™ -2] 
s 


(25) 


dE^(PL)  =  ndP^-b’)!^^  MdP/t-)(S"-T^"F  cos^e’)] 


jn  T  r-r>C^Q 


e-jks 

s 


+  0[m"^,m-3] 


(26) 
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(b)  cl'P'g(Q' )  case: 


dE"(P,  )  =  dP-(Q')sin0’[H*+T  cose’]  +  0[m"^] 

e  L  4  IT  ^  ^  s 


(27) 


dE^(P,  )  =  dP„(Q')sin0’T  F  +  0[m-2] 

e  L  4  IT  s  o  <;  *• 


(28) 


where 


=  g(5i)e  ^ 


S  =  — g(^  )e 
mt(Q  )  i 


(29) 

(30) 


and 


Ki  =  -m)i(Q'  )cos0'' 


(31) 


» (Q') .  — ;W'l 

^  (l+T^2cos2gTp73 


(32) 


The  angle  0’  is  defined  by  n'*s  =  cose’  as  shown  in  Figure  5,  Also, 


F  =  S^-H^cose’ 
I+T^^cos^i 


(33) 


as  defined  in  Table  II. 
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TABLE  II 


(FOR  LIT  REGION) 


KOT  M  «;,e«IC 


F 


D.  Pattern  Factors 

The  solutions  for  short  magnetic  or  electric  dipoles  have  been 
given  in  part  (B)  anu  (C).  One  approach  to  analyze  an  extended  aperture 
or  linear  antenna  problem  is  to  integrate  the  above  solutions  over  the 
source  distribution,  if  it  is  known.  This  is  an  application  of  the 
superposition  theorem,  and  one  approximates  the  source  distribution  by 
an  array  of  short  magnetic  (or  electric)  dipoles  on  the,  conducting 
surface.  This  is  an  accurate  solution,  however,  rather  tedious.  A  more 
efficient  approach  is  to  modify  (j^(Q')  as  shown  in  Reference  [6]  such 
that 

(a)  in  the  shadow  region: 

r  kB  «  "IP  kA  «  • 

p  ^  ^  ZB  cos(7-(Pni.t'))  sin(^Pn,.b') 

"Iff  ::  ?  kA  - 

Ll-(~(Pm-t'))  J  L  T  Pm-b*  J 

Pe  =  n'[l-cos(kL)]  (35) 
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(b)  in  the  lit  region; 


P  2B 

kB  j  ' 

cos(7^in9  (Pn,.t')) 

r  kA  .  -  -  1 

sin(7^ine  Pn,.b') 

ni  tn  ^ 

IcB  ? 

l-{— sine  (Pm-f)) 

kA  i  *  » 

fsine  (Pn,.b') 

(36) 


Pe  =  n'  cos(kL'ri' »s)-cos(kL) 


l-(n' 


(37) 


Here  P^,  =  unit  vector  in  the  direction  of  magnetic  current  moment, 

A,B  =  the  length  of  the  short  and  long  sides  of  the  slot,  and 
L  =  the  length  of  the  monopole. 

It  is  noted  that  L  is  not  to  exceed  a  quarter  wavelength  for  the 
solutions  to  be  /alid. 


E.  Ellipsoid  Surface  Parameters 

The  formulations  in  section  (8)-(D)  are  used  to  solve  for  the 
radiated  fields  of  antennas  mounted  on  an  ellipsoid.  Using  the 
ellipsoid  geometry  shown  in  Figure  6,  the  surface  is  defined  by 

R(9,(>)  =  R(9,4i)  sin9  cosiji  x  +  R(9,4>)  sine  sin()>  y 

A 

+  R(9,*)  cose  z  (38) 

or, 

4-  AAA 

R(VetVr)  “  ^  cosVg  cosVp  x  +  b  cosVg  sinVp  y  +  c  sinVg  z  (39) 
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^dnV  *  ^  cos  6  dnd 

®  c  sin  0 


tanV_  =  ?  sin 


cos  4. 


The  Vp  and  Vg  parameters  are  introduced  because  of  those  convenience  in 
analyzing  elliptic  geometries. 

Considering  a  ray  which  propagates  along  a  geodesic  path  O'Q  on  the 
ellipsoid  surface  as  shown  in  Figure  7,  the  three  unit  vectors  t,  n,  and 

A 

b  are,  as  defined  earlier,  the  geodesic  tangent,  outward  surface  normal 
and  bi normal  at  any  point  along  the  geodesic  path.  The  outward  surface 
unit  normal  (n)  is  obtained  from 


n  =  ''r 


’^Vp  ^Vg 


where 


D  =  3R 
Ve  Wp 


_  -a  sin  Vg  cos  Vp  x  -  b  sin  Vg  sin  Vp  y  +  c  cos  Vg  z 


_3R^  =  -a  cos  Vg  sin  Vp  X  +  b  cos  Vg  cos  Vp  y 


Then, 


be  cos^Vg 
[b^c^  cos^V 

e 


*  2  “ 
cosVp  X  +  ac  cos  Vg  sinVp  y 

cos^v^  +  a^c^  cos^V-  sin^V^ 
r  e  r 


+  ab  cosVg  sinVg  z 

+  a^b^  cos^v  sin^V  ]  V? 
e  e-^ 


cosVe  cosVp  X  +  cosVg  sinVp  y  + 

a  _  b 

A 


A 

sinVg  2 
c 


(42) 


where  A  = 


The  normal  curvatures  on  the  surface  are  evaluated  by  introducing  the 
first  and  second  fundamental  forms  of  differential  geometry  [5]  such 
that 


„  _  L  dVp  +  2M  dVpdVg  +  NdVg 

n  ■  - r - j 

E  dVp  +  2F  dVpdVg  +  GdVg 


(44) 


where 


A 

“  SrVr 

•  N  , 

M  =  R, 

+ 

-► 

E  3  R 

C  «yp 

^Vr  ’ 

F  =  R, 

1  =  -  n. 

VrVe  ’  ^  ^  ~  '^VeVe  *  ^ 


*  ■¥ 

m  P 

Vr  '^Ve 


^e  •  ^e 
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It  can  be  shown  that 


Ru  v  "  -a  cos  V  cos  V  x  -  b  cos  V  sinV  y-csinV  z 

e  e  e  r  e  r  e 

Rw  w  =  a  sin  V  sin  V  x  -  b  sin  V  cos  V  y 
e  r  e  r  e  r 

>  ^  A 

Ry  V  =  -a  cos  V  cos  V  x  -  b  cos  V  sin  V  y  . 

r  r  e  r  e  r 

After  some  algebraic  manipulation,  one  obtains 

2 

,  ^  cos  Vg 
~  ~A 

M  =  0 

N  =  i 

A 

E  =  a2  cos^  Vg  sin2  Vp  +  b^  cos^  Vg  cos^  Vp 
F  =  (a^  -  b2)  sin2  gin^  Vp 

G  =  a^  sin2  Vg  cos^  Vp  +  b^  sin^  Vg  sin2  Vp  +  c^  cos^  Vg 

A  pair  of  orthogonal  directions  exists  for  which  curvature,  K, 
assumes  maximum  and  minimum  values,  i.e.,  principal  directions 

A  A 

represented  by  two  unit  vectors  and  rp*  Two  extreme  values  of  K 
corresponding  to  the  above  directions  are  called  principal  curvatures 
denoted  by  Kj  and  K2. 


?A 


Mean  curvature:  Kik|  -  *^1  ^  ^2  =  EN  -2MF  +  LG 

2  2(EG-F^) 

Gaussian  curvature:  Kg  =  Kj  <2  =  LN-M^ 

EG-F^ 

Thus,  the  principal  curvatures  and  K2  are  given  by 

Kl.  2  =  Km  ±  • 

A  A 

The  two  principal  directions  (ti,t2)  are  given  by 
t1  =  1  fl  Rv  +  a  Rv  1 

n  L 

,2  =  ^  [s  R*v^  .  1  Rv  1 


(45) 


(46) 


where 


a  =  L  -  KiE  0  =  M  -  K2F 

KiF  -  M  ,  K2E  -  L 


Yi  *  (E  +  2aF  +  a^G) 


Y2 


=  (e^E  +  28F  +  G) 


(47) 
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However,  it  is  noticed  that  ~  and  K2  ~  within  the  significant 
energy  region  for  most  practical  cases.  That  indicates  approximate 
values  of  the  two  principal  curvatures,  ~  and  K2  ~  ^/q,  are  good 
enough  to  be  used  for  the  geodesics  on  the  ellipsoid  surface  for  our 
radiation  consideration. 

Thus 


K 


1 


L 

E 


cos^  Ve 

A  [a^  cos2  sin2  y  +  cos^  cos^  v^] 
G  r  6  r 

=  (A[a2  sin2  +  b^  cos^  Vp]) 


(48) 


K 


2 


■»  1 
G 


_  1  _ _ _ 

A  [a2  sin2  y  cos^  V  +  b^  sin^  sin^  +  c^  cos^  V^]  • 

e  r  e  r  e 


(49) 


It  is  noticed  that  r  =  1_  and  r  =  as  found  in  Table  I.  For  most 

1  Ki  2  K2 

practical  cases,  a  and  b  within  the  significant  energy  region  become 

very  small.  Thus,  one  may  use  Ry  and  Ry  as  the  principal  surface 

r  e 
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directions,  which  are  expressed  by  tj,  respectively. 

A  A  A  A 

If  6'  denotes  the  angle  between  t  and  ti,  then  t  =  ti  cos  0*  +  t2  sin  0’. 
From  Euler's  theorem,  the  normal  curvature  along  the  geodesic  path  is 
specified  by 


Kg  =  Kj  cos2  0'  +  K2  sin2  0'  (50) 

with  the  radius  of  curvature,  pg,  being  Vkg» 

The  torsion  term  (Tq)  introduced  in  Section  (B)  is  given  by 
To  =  T  •  Pg  (51) 

where  the  surface  torsion  is  given  by 


T  =  sln|  0J.  (K^  .  Kj) 


(52) 


with  K\  and  <2  being  defined  in  Equations  (48)  and  (49). 


CHAPTER  III 

NUMERICAL  TECHNIQUE  AND  PERTURBATION  METHOD 
A.  Introduction 

It  is  seen  that,  for  an  antenna  mounted  on  an  ellipsoid,  the 
geodesic  paths  associated  with  the  GTD  solution  in  the  shadow  region  are 
extremely  complex.  An  elaborate  method  employed  calculus  of  variations 
to  calculate  the  geodesic  paths  which  resulted  in  a  very  Inefficient 
solution  is  presented  in  Chapter  IV.  An  efficient  numerical  approach  is 
examined  in  this  chapter  with  the  ellipsoid  simulated  by  an  elliptic  cone 
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or  elliptic  cylinder  model.  Since  the  elliptic  cone  and  elliptic 
cylinder  are  developed  surfaces,  geodesics  are  efficient  to  solve. 


B.  Surface  Geodesics 

The  geodesics  on  the  elliptic  cylinder,  elliptic  cone  and  ellipsoid 
are  presented  in  this  section. 

a)  Elliptic  Cylinder  Case: 

The  elliptic  cylinder  geometry  used  for  this  study  is  shown  in 
Figure  8(a).  Since  the  elliptic  cylinder  is  a  developed  surface,  the 
geodesic  path  Q'Q  is  a  straight  line  on  the  unfolded  planar  surface.  As 
shown  in  Figure  8(b),  the  geodesic  unit  tangent  t  is  given  by 


t  =  Vp  cos  Y  +  z  sin  Y 


(53) 


A  A  A 

where  \j  =  -a  sin  Vp  x  +  b  cos  Vp  y 

^  (a^  sin^  +  b^  cos^  V^)l/2 


(54) 


It  is  noticed  that  Vp  and  z  are  the  two  principal  directions  on  the 
elliptic  cylinder  surface.  For  a  given  geodesic  Q'Q,  one  can  see  that  y 
is  a  constant  along  the  geodesic  path, 
b)  Elliptic  Cone  Case: 

Consider  a  ray  propagates  along  a  geodesic  path  Q'Q  on  the  elliptic 
cone  surface  as  shown  in  Figure  9(a).  It  is  a  straight  line  on  the 
unfolded  planar  surface  as  shown  in  Figure  9(b).  It  is  noticed  that  Vp 

A 

and  tg  are  the  two  principal  directions  on  the  surface.  The  geodesic 
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r 


unit  tangent  t  is,  then,  represented  by 

A  A  A 

t  =  tp  COS  6  +  tg  sin  B 

where 


tg  =  -Xg  sin  5  +  z  cos  5 

Xg  _  X  a '  cos  Vp  +  y  b'  sin  Vp 
a'2  cos^  Vp  +  b'2  sin^  Vp 

and  6  is  the  half  cone  angle  as  shown  in  Figure  9(a).  Note  that 
longer  a  constant  along  the  geodesic  path  Q'Q. 

(c)  El lipsoid  Case: 

Recall  that  the  ellipsoid  surface  is  defined  by 

R(Ve>'^r)  "  ^  '^e  Vp  x  +  b  cos  Vg  sin  Vp  y  +  c  sin  Vg 

with 

tanV  =  ^  9 

®  c  sin  9 

and 

tan  V  =  ^  sin  .)> 
b  cos  t 

From  Section  II-(E),  it  was  found  that  the  principal  directions 
*► 

given  by  R^g  and  Ryp.  As  shown  in  Figures  10(a)  and  (b)  the  uni 
in  the  principal  directions  can  be  represented  by 


(55) 

(56) 

(!>7) 

B  is  no 


were 

t  vectors 
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X 


SOURCE 


e  10(a).  Geodesic  path  on  an  ellipsoid. 
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Figure  10(b).  Elliptic  profile. 
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tp  ~  Rv  /  iRv  1 

r  r 


— a  cos  Vg  sin  x  +  b  cos  Vg  cos  Vp  y 
cos^  Vg  sin2  Vp  +  cos^  Vg  cos^  Vp 


tg  ~  Rw  /  \R^J 
e  e 


(58) 


~  -a  sin  Vg  cos  Vp  x  _  b  sin  Vg  sin  Vp  y  +  c  cos  Vg  z 
sin^  Vg  cos2  Vp  +  b2  sin^  Vg  sin2  cos^  Vg 


(59) 


If  0  denotes  the  angle  between  tp  and  the  geodesic  unit  tangent  t,  then 

A  A  A 

t  =  tp  COS  6  +  tg  sin  6  .  (60) 

which  is  identical  to  the  form  used  for  the  elliptic  cylinder  and 
elliptic  cone  geodesics.  This  suggests  that  one  might  be  able  to  develop 
a  perturbation  solution  which  gives  a  simplified  form  for  s  on  an 
ellipsoid  using  the  elliptic  cylinder  or  elliptic  cone  expressions  for  0. 


C.  Elliptic  Cylinder  Perturbation 

In  order  to  solve  for  the  geodesics  on  an  ellipsoid,  the  elliptic 
cylinder  perturbation  technique  is  used  when  the  source  is  located  at 
0s  =  90®.  Note  that  the  source  position  is  assumed  to  be  in  the  ((>$  =  0 
plane. 
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Recalling  that  y  is  a  constant  along  a  given  geodesic  path  on  an 


elliptic  cylinder,  one  obtains  a  geodesic  equation  given  by 


tan  Y  =  •  (61) 

The  elliptic  cross-section  in  the  Xg  -  z  plane  is  described  by 


The  elliptic  cross-section  in  the  x’-y'  plane  is  described  by 

x'  =  a  cos  Vg  cos  Vp 
y'  =  b  cos  Vg  sin  Vp 
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therefore. 


At  the  diffraction  point  Q,  the  radiation  direction  (fi^,  should 
coincide  with  the  geodesic  tangent  t  given  in  Equation  (60).  Thus, 


t  =  X  sin  9t  cos  +  y  sin  9t  sin  ^  *  z  cos  9^ 

A  A 

=  tp  cos  Y  +  tg  Sin  y  (65) 

In  order  to  trace  the  geodesic  path  from  the  source  location  to  the 
diffraction  point,  one  can  always  assume  a  diffraction  point  at  the 
source  (Vg,  Vp)  =  (0,0)  with  the  radiation  direction  (e  ,  ^  ^)  and 
gradually  add  increments  {n&i,  A(t>t)  until  the  desired  radiation  direction 
(9t*  't*t)  reached. 
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After  the  geodesic  path  is  determined,  it  is  necessary  to  check  the 


significant  energy  region  by  calculating  the  Fock  parameter.  The  Fock 
parameter  C  of  Equation  (12)  is  given  by 

5  =  /  1.  1/3 

Q'  pg  ♦  (*56) 

where  the  integral  is  evaluated  along  the  geodesic  path.  Note  that  t  is 
the  arc  length  along  the  geodesic  and  given  by  either  of  the  following 
equations; 


Therefore 
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or 

Vr  ,  1 

=  1  I  1_  3  . 

cos  y  0  pg  \~2~ J  (f’S) 

\/  cos^  sin2  v'  +  cos^  v  cos^  dv'  • 

V  e  r  e  r  r 

where  pg  =  1  /  (kj^  cos^  y  +  *<2  sin^  y)  and  kj,  k2  are  given  in  Equations 
(48)  and  (49). 

D.  Elliptic  Cone  Perturbation 

When  the  source  is  not  located  at  the  mid-section  (65  *  90°),  the 
elliptic  cone  perturbation  method  is  used.  The  geodesic  path  on  an 
ellipsoid  using  the  elliptic  cone  perturbation  method  is  shown  in  Figure 
11  and  the  associated  unfolded  surface  is  shown  in  Figure  12.  If  y  and  s 

A  A 

denote  the  angle  between  t  and  tp  at  O'  and  Q  respectively,  it  is  seen  in 
Section  III-(B)  that  they  are  not  the  same  as  in  the  elliptic  cylinder 
case. 

In  fact, 

6  =  y  -  a  .  (89) 

The  calculation  of  a  will  be  discussed  in  detail  in  this  section. 

^  A  A  A 

r  =  X  a  cos  cos  Vp5  +  y  b  cos  sin  +  z  C  sin  Vg^ 

AAA 

♦  + 

t  =  3r  =  -X  a  sin  Vgs  cos  Vpj  -  y  b  sin  Vgj  sin  Vp5  +  z  c  cos  Vgg 
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Geodesic  path  on  the  perturbed  elliptic  cone. 
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1 


therefore,  one  can  show  that 


x-a  cos  Vgs  cos  Vps  _  y-b  cos  Vgj  sin  Vpj 


-a  sin  Vg5  cos  Vpj 


-b  sin 


'es  ^rs 


z-c  sin  Vgs 
c  cos'  Vgs 


For  x=y=0,  z  =  c  esc  Vgj 


Accordingly, 


rs  =  [a2  cos2  Vgs  cos^  Vps  +  b2  cos^  Vgs  sin?  Vpg 
+  (c  sin  Vge  -  c  CSC  Vgs)2] 


=  [a?  cos?  Vps  +  b?  sin?  Vps  +  cot^  Vgj] 


where  Zj  =  c  sin  Vgj 
as  =  a  cos  Vgj 


bj  =  b  cos  Vgs 


From  Figure  13(b)  and  (c). 


1  =  X  as  cos  Vps  +  y  bs  sin  Vps. 


at 


-  -X  as  sin  Vps  +  y  bs  cos  Vps 


From  the  previous  calculation,  one  can  define  a  new  vector  U  as 


following: 

IJ  = 


(^)  '  (^) 
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Figure  13(a).  Elliptic  profile 
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+  Z  C  cos  Vgs) 

^  ^  ^  /s 

^  (-X  85  sin  Vps  +  y  bg  cos  Vrs),^  (-X  a  sin  Vgs  cos  Vps  -  y  b  sin  Ves  sin  V^g 
[a^  sin2  cos^  [a?  sin^  cos  TIT— rsr  sin^  Vg^  sin^ 

+  c^  cos^  V  ]  1/2 
es 


A  A 

=  X  bg  c  cos  Ves  cos  ^rs  +  y  ^s  c  cos  Vgs  sin  V^s  +  z  ag  b  sin  Ves 

__ 

[a2  sin2  V  +  b^  cos^  V  ]  [a^  sin^  v„  cos^  V_.  +  b^  sin^  \i  sin^ 

^  I  j  o  » b  “s  •  s  c5  rs 

V2 

+  c^  cos^  V^^] 
es 


The  magnitude  of  vector  U  becomes 


sin  0 

R 


=  sin  9 

R 


From  Figure  13(c),  one  obtains  the  following  equation; 
da  = 

f's 

Therefore  a  is  given  by 


a 


opp  ??  1/2 

cos  Ves  [a  b  sin'  Ves  +  c '  cos  Ves  (a  sin  Vp  +  b  cos  Vp)] 

c  [c2  cos2  V  +  sin2  y  (a^  cos^  v'  +  b^  sin^  V')]  17? - "r 

s  es  es  r  r 


(70) 
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r 

! 


i 

! 


where  r 

s 


=  [a^  cos^  V'  +  sin^  V'  +  zf  cot^ 
s  r  s  r  s 


''es^l 
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To  obtain  a  diffraction  point  (Vg,  V^)  for  a  radiation  direction  (e^, 
one  can  always  assume  the  first  diffraction  point  is  at  the  source 
(Ve»  Vp)  =  (Vj,  0)  with  the  radiation  direction^’',  +.  ^^and  gradually  add 
the  increments  (A8t,  until  the  final  radiation  direction  (0^,  4)^) 
is  reached. 

The  geodesic  arc  length  is  obtained  from  either  of  the  following 
equations ; 

t  cos  Y  =  (rg  -  Sg)  sin  a  ,  or  (71) 

l  sin  Y  =  r  -  (rg  -  Sg)  cos  a  .  (72) 

Therefore 

di  =  r  cos  Y  da  »  0^ 
cos2(Y-a) 


dt  =  - 


1 

sTrTTr^ 


dV7 


('’c 


Se) 


dV, 


The  Fock  parameter  ^  is  obtained  by  integrating  along  Vp  or  Vg,  i.e.. 


?  =  r  cos  Y 


1 

COs2(Y-a' ) 


da' 


=  r  cos 


1 

COS^( Y-a) 


da  dVp 

3V7 


(73) 
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where 


2w2  ...2  „  .  .2  „.2  .  ,2  .„2  „  j^/2 


da  -  cos  Ves  [a  ^  ^es  ■*■  ^  ^es  '^rs  '*■  ^  '^rs 

- 172 - 


rs 


r^  [c2  cos2  y  +  sin^  y  (a^  cos2  y  +  sin2  v„)] 
s  es  es  '  rs  rs ' 


or 


C  =  / 

0  Pg 


(>) 


I 

3  1  d  (r  .  ^  ) 

sTnrP^  W;  ^  s 


dy. 


(74) 


where  = 


1 


and 


Ki  cos2  Y  +  K,  sin2  y 


\  ^  WV/-9  I  »  IX^ 

Kj,  <2  are  given  in  Equations  (48)  and  (49). 


CHAPTER  IV 

EXACT  GEODESIC  PATHS  OF  AN  ELLIPSOID 


More  elaborate  numerical  method  employing  calculus  of  variations  to 
calculate  the  geodesic  paths  on  an  ellipsoid  is  studied  in  this  chapter 
in  order  to  check  the  validity  of  the  previous  perturbation  solutions. 
Using  rectangular  coordinates,  an  ellipsoid  can  be  described  by 

4  4  4  =  1 

a?  ^  c^ 

where,  without  loss  of  generality,  it  is  assumed  that  c  >  b  >  a  >  0. 

The  key  to  the  derivation  of  the  geodesic  path  solution  of  an  ellipsoid 
is  to  find  a  coordinate  system  which  is  orthogonal  on  the  ellipsoidal 
surface. 
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Consider  the  following  three  equations: 


lO  TO 


1  .  a2  >  C 

1  ,  b2  >  Ti  >  a2 
1  .  c?  >  T  >  b2 


(76) 

(77) 

(78) 


which  are  respectively  of  an  ellipsoid,  a  hyperboloid  of  one  sheet  and  a 
hyperboloid  of  two  sheets,  all  confocal  with  the  ellipsoid  of  Equation 
(76).  The  variables  u^  =»  c.  u^  =  n,  u^  =  t  are  called  ellipsoidal 
coordinates.  The  transformation  to  the  rectangular  coordinates  is 
obtained  by  solving  Equations  (76),  (77),  and  (78) 


simultaneously  for  x,  y  and  z,  such  that 


(79a) 

(79b) 

(79c) 


In  terms  of  the  ellipsoidal  coordinates,  the  displacement  vector  dr 


can  be  written  as 


dr  =  du^  (j(j2  +  3r  (jy3 


3uJ 


3u^ 


n  *►  45  "►  O 

=  a^du^  +  a2du^  +  a3du-5 

Then  the  length  of  a  line  element,  denoted  by  ds,  is 


(80) 


-  +  +  3  3  4.  +  ,  . 

{ds)2  =  dr  •  dr  =  y  ai  •  du^  duJ 
i=l  j=l 


+•  + 


where  q. .  =  a.  •  a . 

1  J 


3  3 

I  I  9ij  duJ 
i=l  j=l 


3x  3x  +  3y  3y  +  3z  3^ 

Su"*  3u3  3u^  3uJ  3u'  3uJ 


(81) 


(82) 


It  can  be  shown,  using  Equations  (79a),  (79b)  and  (79c),  that 


9ij  =  0» 

if  i*j 

(83a) 

911  =  i 

(T-0(n-0  . 

(83b) 

T 

(a2.0(b2-C)(c2-c) 

922  =  1 

(C-r))(T-n) 

(83c) 

T 

(a2-n)(b^-n)(c2-r\) 

933  =  1 

(n-T)(c-T) 

(83d) 

T 

(a^-x) (b^-t) (c^-t) 

By  substituting  Equations  (83a)  through  {83d)  into  Equation  (81),  one 
obtains 
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(<is)2  .  ^ 


. —JLlliKvSL _  (()?)2  +  ■ 


^  ..(.^-lIU---^l -  {dx)A 

(a2-T)(b2-T)(c2-T)  ; 


(84) 


For  the  ellipsoid  of  Equation  (75),  c=dc=0.  Then  Equation  (84)  becomes 


(ds)2  =  1.  (T-n) 


n(dn)2 

( n-a2 ) ( n-b^ ) ( ) 


jlM 

iTt^ 


(T-a2)(T-b2)(T-c2) 


(85) 


or 


nn  ^ 

( n-a  2 )  ( n-  b^yi  ri-c^ ) 


(T-a2)(T-b2)(T-c2)f 


(86) 


where  „•  -  dn  and  dr  is  assumed  to  be  positive.  Let  h(n,n’,T)  denote 

17 


the  integrand  of  Equation  (86),  i.e.. 


V2 


h  s  (r-n) 


1/2 


nn 


<2 


( n-a2)  ( n-b2)  (ri-c^)  ( t-a  2 )  ( r-b^ )  ( V-c^ ) 


(87) 


Then  using  the  calculus  of  variation  technique,  it  can  be  shown  that 
the  geodesic  path  satisfies 
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t 


The  Equation  (88)  becomes  a  complete  differential.  The  derivation  of 
the  geodesic  path  solution  is  straight  forward  from  Equation  (91). 

It  is  obvious  that 

n'(T-n)  i  +  n  -  b^  =  -8  (92) 

n  3n' 

where  0  is  a  constant. 

By  combining  Equation  (87)  with  Equation  (92),  one  obtains 

-n(dn)^ _  =  _ -T(dT)^ _ 

( n-a^ ) ( n-b^ ) ( n-c^ ) ( n-b2+ e)  ( r-a^) ( T-b^ ) ( t-c^) ( t-b^+e) 

(93) 

Since  c2>T>b2>n>a2,  it  is  obvious  that  T-b2+8>0>n-b2+8  or  b2-n>8>b2-T. 
Equation  (93)  is  the  geodesic  path  solution.  However,  it  is  more 
convenient  to  make  the  following  changes  of  variables. 

n  =  a2  sin2  (|)  +  b^  cos2  ^  ,  and  (94a) 

T  =  b2  cos2  ij;  +  c2  sin2  ifi  ,  (94b) 

In  terms  of  ^  and  iji,  the  geodesic  path  solution  in  Equation  (93) 


can  be  rewritten  as 


and  in  integral  form  as 


V2 


[b2-_jb£-a£)_sin£_^ 


[(c^-b^)  +  (b^-a^)  sin^  +  (b^-a^)  sin^  iji] 


1/2 


/  1  _  b^  +  (c^-b^)  sin^  ji  _ 

'i'l  I  [(b2-a2)  +  (c2-b2)  sin2  ^][0  +  (c^-b^)  sin^  ’J'j 


d4i 


dt)^ 


(96) 


where  b2-a2  >  g  >  b^-c^.  Note  that  the  absolute  values  of  dit*  and  d\|;  are 
used  in  Equation  (96). 

To  define  the  ranges  of  it*  and  tp,  it  is  necessary  to  go  back  to 
Equations  (79a),  (79b)  and  (79c).  With  5=0  and  in  terms  of  41  and 
they  can  be  rewritten  as 


V2 

I _  cos  <|)  (b2  cos^  ii/  +  sin  2  ^  .  a2)  (97a) 


T/2 


(c2.a2) 


y  =  b  sin  <t»  sin  41 


(97b) 


V2 

z  =  ^  cos  4»  (c2-a2  sin2  4,  -  b2  cos2  4,)  .  (97c) 

\  ^^2 
(c2-a2) 


If  the  geodesic  path.  Equation  (96),  crosses  the  curve  4)=0  or  tt,  then 
e<0  and  if  it  crosses  the  curve  \|;=0  or  n,  then  R>0.  Thus  to  ensure  the 
continuity  of 4)  and  4;  along  the  geodesic  path,  the  ranges  of  41  and  4»  will 
be  defined  as  follows: 
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TT  for  B  <  0 

and  -IT  ^  ij)  £  2tt,  0  £  41  ^  n  for  b  >  0. 

Figure  14  illustrates  the  4>  and  4;  curves  as  projected  onto  the  x-z 
plane.  Equation  (96)  is  the  geodesic  path  solution  of  an  ellipsoid 
which  is  employed  for  obtaining  more  elaborate  geodesic  paths.  When  the 
geodesic  path  includes  one  of  those  four  points,  Pj,  P2,  P3  and  P4  in 
Figure  14,  then  6  must  equal  zero  and  the  integrals  in  Equation  (96) 
diverge  at  those  points.  Thus  if  (4>2.  '^2)  =  Pi .  i  =  1,  2,  3,  4,  then  it 
is  necessary  to  replace  (4)2,  412)  ('^*2  'l'2  where  A41  -  0  and 

A4»  -  0.  The  geodesic  path  between  the  two  points  (4.1,  411)  and  (4)2,  4<2) 
can  be  determined  from  Equation  (96)  by  first  determining  the  value  of 
B.  Since  there  are  absolute-value  signs  attached  to  the  differential  d4) 
and  d4i  in  the  geodesic  path  solution,  it  is  important  to  know  how  the 
variables  41  and  4;  change  from  4>i  to  412  and  4;^  to  412,  respectively. 

CHAPTER  V 
RESULTS 

Geodesics  for  a  source  mounted  on  an  ellipsoid  can  b^ .analyzed 
precisely  by  computing  the  geodesic  path  defined  by  the  surface 
parameters  (0q,  4)g)  and  the  geodesic  tangent  defined  by  the  radial 

A  A 

vector  direction  (0t,  4>t)  as  shown  in  Figures  15(a)  and  (b).  The 
geodesic  path  indicates  the  actual  diffraction  point  location  on  an 
ellipsoid;  whereas,  the  geodesic  tangent  indicates  the  radiation 
direction  at  the  corresponding  diffraction  point. 

To  show  the  validity  of  the  elliptic  cylinder  perturbation 
solution,  the  source  is  placed  at  05=90°  and  the  geodesic  paths  and 
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\I/B  ir 


-  <#>=  CONSTANT 

- ^  =  CONSTANT 


Projection  of  and  ip-  curves  onto  the  xz  plane 
of  an  ellipsoid. 


SOURCE 

LOCATION 


geodesic  tangents  associated  with  this  source  location  are  calculated  as 
shown  in  Figures  16(a)  through  (d)  and  17(a)  through  (d),  respectively. 

In  each  figure,  the  elliptic  cylinder  perturbation  solutions  are 
compared  with  the  exact  solutions.  Looking  through  those  figures,  one 
can  see  that  the  geodesic  paths  and  the  geodesic  tangents  of  both 
methods  coincide  with  each  other  within  the  significant  energy  region 
close  to  the  source.  This  coincidence  in  the  significant  region  can  be 
checked  more  precisely  by  calculating  the  Fock  parameter  (^)  along  each 
geodesic  path  as  shown  in  Figures  16(d)  and  17(d).  Actually,  the  Fock 
functions  associated  with  the  solutions  drop  more  than  20db  as  the  Fock 
parameter  (c)  reaches  2.5  in  the  deep  shadow  region.  This  clearly  shows 
the  significant  portion  of  the  surface  as  discussed  in  Chapter  I.  It  is 
also  noted  that  the  geodesic  paths  on  the  2xx3xx20x  ellipsoid  shown  in 
Figures  16(b)  and  17(b)  show  better  agreement  than  those  on  the  2xx3x4x 
ellipsoid  shown  in  Figures  16(a)  and  17(a).  Thus,  the  solution  becomes 
more  correct  for  shapes  which  are  more  cylindrical  .  So  one  would  expect 
that  this  approximate  solution  will  work  quite  successfully  for  more 
realistic  aircraft  and  missile  shapes  which  tend  to  be  more  cylindrical. 

The  elliptic  cone  perturbation  solutions  can  be  examined  by  placing 
the  source  at  05=30®  as  shown  in  Figures  18(a)  through  (c),  where  the 
geodesic  paths  of  the  elliptic  cone  perturbation  solutions  are  compared 
with  those  of  the  exact  solutions  for  three  different  ellipsoid  geomet¬ 
ries.  One  should  note  the  good  agreement  between  the  two  results.  The 
significant  energy  region  of  the  geodesic  paths  is  shown  using  the 
calculated  Fock  parameters  (^)  in  Figure  18(c).  It  is  also  noted  that 
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-  PERTURBATION 

-  EXACT 

r  =  280‘’ 


Figure  16(a).  Geodesic  paths  defined  by  the  surface  parameters 
(0Q»  4^)  ^0*^  3  source  mounted  at  05=90°  on  a  2xx 
3Xx4x  ellipsoid.  Note  that  y  is  the  angle  between 
the  geodesic  tangent  t  and  one  principal  direction 

A 

tp  at  the  source  location. 
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-  PERTURBATION 

-  EXACT 


Figure  16(b).  Geodesic  paths  defined  by  the  surface  parameters 
(9q*  a  source  mounted  at  0,^=9ri'’  on  a  2xx 

3Xx?0X  el  1 ipsoid. 
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Figure  17(a).  Geodesic  tangents  defined  by  the  radial  vector 

direction  (9^,  <j>t)  for  a  source  mounted  at  85=90° 
on  a  2xx3xx4x  ellipsoid.  Note  that  y  is  the  angle 
between  the  geodesic  tangent  t  and  one  principal 
direction  tp  at  the  source  location. 
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-  PERTURBATION 

- -  EXACT 


Figure  17(b).  Geodesic  tangents  defined  by  the  radial  vector 

direction  (9^,  (|>t)  for  a  source  mounted  at  85=90° 
on  a  2xx3Xx20x  ellipsoid. 
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<p^  {  DEGREES) 

Geodesic  tangents  defined  by  the  radial  vector 
direction  (9t>  for  a  source  mounted  at  65=90 
on  a  2Xx4xx20X  ellipsoid. 


-  PERTURBATION 

-  EXACT 


o'  •  I  •  ■  I  ■  •  I 

90.  120.  150.  180. 

(  DEGREES) 


Figure  17(d).  Geodesic  tangents  defined  by  the  radial  vector 

direction  (9^,  (t^^)  for  a  source  mounted  at  85=90 
on  a  4Xx6\x40X  ellipsoid. 


-  PERTURBATION 

-  EXACT 


Figure  18(b).  Geodesic  paths  defined  by  the  surface  parameters 


(®Q»  ♦q)  ^or  a  source  mounted  at  05=30°  on  a  ?.\k 
3Xx8x  ellipsoid. 


PERTURBATION 


EXACT 


</)q  (  DEGREES) 

Figure  18(c).  Geodesic  paths  defined  by  the  surface  parameters 
{9q»  i>())  3  source  mounted  at  65=30°  on  a  4xx 

6xx40x  el  1 ipsoid. 
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-  PERTURBATION 

- EXACT 


Figure  18(d).  Geodesic  paths  defined  by  the  surface  parameters 


(9g,  (frg)  for  a  source  mounted  at  05=86°  on  a  4xx 
6Xx40x  el  1 ipsoid. 


the  geometry  of  the  ellipsoid  used  in  Figure  18(c)  is  more  realistic  in 
terms  of  simulating  an  aircraft  or  missile  fuselage  than  those  used  in 
Figure  18(a)  and  (b).  To  make  sure  that  the  solution  switches  correctly 
between  the  elliptic  cylinder  and  the  elliptic  cone  perturbation  method, 
the  source  is  placed  at  95=86°  and  the  geodesic  paths  are  calculated  as 
shown  in  Figure  18(d).  Comparing  Figure  18(d)  with  Figure  16(d),  one 
can  see  that  the  geodesic  paths  of  the  elliptic  cone  perturbation 
solution  are  very  close  to  those  of  the  elliptic  cylinder  perturbation 
solution  for  that  source  location. 

In  Figures  19(a)  through  (p),  the  geodesic  tangents  obtained  by  the 
elliptic  cone  perturbation  solution  are  compared  with  those  of  the  exact 
solutions  for  several  different  source  locations  and  ellipsoid 
geometries.  As  one  moves  the  source  location  away  from  the  two  ends  of 
the  ellipsoid,  one  can  see  that  the  discrepancies  between  perturbation 
solutions  and  exact  solutions  become  smaller  as  shown  in  the  sequence  of 
figures,  i.e..  Figures  19(a),  (b),  (c)  and  (d).  In  order  to  provide 
more  information  about  the  applicable  geometries  and  source  locations, 
several  different  geometries  and  source  locations  are  examined  i 
Figures  19(e)  through  (0).  For  the  more  realistic  case,  i.e.,  the 
4Xx6Ax40A  ellipsoid,  with  the  source  location  9s=30°  is  plotted  in 
Figure  19(p).  Although  one  can  see  small  discrepancies  for  the  rays 
toward  the  tips  of  the  ellipsoid,  i.e.,  lines  for  y=80°  and  y=280°  in 
the  Figure  19(p),  they  only  happen  when  the  caustic  effects  come  into 
play.  In  the  caustic  region  where  virtually  an  infinite  set  of  rays 
have  significant  effects  on  the  radiation  pattern,  the  basic  RTD  theory 
fails.  The  study  of  this  caustic  effect  is  beyond  the  scope  of  this 
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study.  If  one  neglects  the  caustic  regions,  the  geodesic  tangents  using 
the  perturbation  model  shown  in  Figure  19(p)  coincide  with  the  exact 
solution  very  well  in  the  significant  region. 


CHAPTER  VI 

SUMMARY  ANH  CONCLUSIONS 


The  object  of  this  study  is  to  develop  an  efficient  numerical 
technique  to  calculate  the  geodesic  paths  of  an  ellipsoid-mounted 
antenna.  The  Geometrical  Theory  of  Diffraction  is  the  basic  approach 
applied  here.  The  curved  surface  diffraction  solutions  are  discussed  in 
Chapter  II,  where  the  creeping  wave  solutions  in  the  shadow  region  are 
of  particular  interest. 

Elliptic  cylinder  and  elliptic  cone  perturbation  methods  are 
presented  in  Chapter  III  to  simulate  the  geodesic  paths  on  an  ellipsoid, 
which  in  turn  can  be  used  to  model  an  aircraft  or  missile  fuselage. 
Because  the  elliptic  cylinder  and  elliptic  cone  are  developed  surfaces, 
the  geodesic  paths  can  be  found  on  an  unfolded  planar  surface. 

In  order  to  show  the  validity  of  the  perturbation  solutions,  more 
elaborate  numerical  method  for  the  geodesic  paths  employing  calculus  of 
variations,  whose  results  are  indicated  as  exact  solutions  here,  is 
also  presented  in  Chapter  IV.  Although  this  method  provides  accurate 
geodesic  paths  on  the  ellipsoid,  it  is  too  complicated  and  inefficient 
to  use  for  practical  radiation  applications.  However,  the  exact 
solution  is  most  appropriate  for  coupling  problems  where  the  exact  path 
is  desired  between  two  known  points  on  the  surface. 
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Figure  19(a).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (e^,  (^)  fur  a  source  mounted 
at  05=30°  on  a  2Xx3Xx4x  ellipsoid.  Note 
that  y  is  the  angle  between  the  geodesic 

A 

tangent  t  and  one  principal  direction  t^,  at 
the  source  location. 
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Figure  19(b).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (0^,  for  a  source  mounted 
at  85=45°  on  a  2Xx3xx4x  ellipsoid. 


72 
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Figure  19(c).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (0-^,  for  a  source  mounted 
at  05=60°  on  a  2Xx3xx4x  ellipsoid. 
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Figure  19(e).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (9-^,  <t>t)  for  a  source  mounted 
at  05=30°  or  a  2xx3xx8x  ellipsoid. 


75 
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Figure  19(f).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (e^,  ij)^)  for  a  source  mounted 
at  05=45°  on  a  2xx3Xx8x  ellipsoid. 


Figure  19(g).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (e^,  for  a  source  mounted 
at  05=60°  on  a  2xx3xx8x  ellipsoid. 
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Figure  19(h).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (9t,  ((►{)  for  a  source  mounted 
at  05=80°  on  a  2xx3xx8x  ellipsoid. 
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Figure  19(k).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (6^,  for  a  source  mounted 
at  05=30°  on  a  2xx3xxl2x  ellipsoid. 
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(  DEGREES) 

Geodesic  tangents  defined  by  the  radial  vector 
direction  (e^,  <t>t)  for  a  source  mounted 
at  05=30“  on  a  2xx3xx20x  ellipsoid. 


Figure  19(m).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (0^,  (}>t)  for  a  source  mounted 
at  9s=45°  on  a  2Xx3xx20X  ellipsoid. 
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Figure  19(n).  Geodesic  tangents  defined  by  the  radial  vector 
direction  <|>t)  for  a  source  mounted 

at  05=60“  on  a  2xx3xx20x  ellipsoid. 
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<p^  (  DEGREES) 

Figure  19(o).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (0t,  4it)  ®  source  mounted 
at  05=80*’  on  a  2xx3ax20X  ellipsoid. 


Figure  19(p).  Geodesic  tangents  defined  by  the  radial  vector 
direction  (9t,  <|>t)  for  a  source  mounted 
at  05=30“  on  a  4xx6Xx40x  ellipsoid. 
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The  geodesic  paths  of  the  perturbation  solutions  are  compared  with 
those  of  the  exact  solutions  for  several  different  geometries  and  source 
locations  in  Chapter  V.  The  comparison  of  both  results  illustrates  that 
the  geodesic  paths  can  be  solved  using  either  numerical  technique; 
however,  the  perturbation  is  much  more  efficient.  In  addition,  one  can 
easily  relate  the  radiation  direction  with  the  desired  geodesic  path 
using  the  perturbation  method.  On  the  other  hand,  one  is  not  sure  which 
geodesic  is  necessary  to  achieve  the  desired  radiation  direction  using 
the  exact  solution. 

This  ellipsoidal  model  will  be  applied  next  to  analyze  antenna 
patterns  for  antennas  mounted  on  aircraft.  The  ellipsoid  will  be  used 
to  simulate  the  fuselage. 
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